



GUIDANCE, NAVIGATION 
AND CONTROL 


U 

< 

to 

i/> 

< 

iS 


MIX CHARLES 

CAMBRIDGE MASSACHUSETTS 02139 


STARK DRAPER 
LABORATORY 




GUIDANCE, NAVIGATION 
AND CONTROL 


Approved: 


Approved: 




fe rt P. F reeman 
Deputy Associate Director 
Charles Stark Draper Lr-Hnratory 




Davifl G . 

Apollo Gfjricjance 
Navigation Prog 



edtor 


Approved: 


pj-jfvL Gc 

RalphLR. Ragan 
Deputy Director 
Charles Stark Draper Laboratory 



E-2648 


CROSS FLOW IN A STARVED 
EHD CONTACT 

by 

Edward P. Kingsbury 
March 1972 


MIX 


CAMBRIDGE, MASSACHUSETTS, 02139 


CHARLES STARK DRAPER 

LABORATORY 


To be presented at the ASLE Annual Meeting 
in Houston, Texas, May 1-4, 1972. 



ACKNOWLEDGMENT 


This paper was prepared under the auspices of DSR Project 
55-23829, which is sponsored by the Manned Spacecraft Center of 
the National Aeronautics and Space Administration through Contract 
NAS 9-4065. 

Publication of this report does not constitute approval by NASA 
of the findings or conclusions contained herein. It is published for 
the exchange and stimulation of ideas. 


© Copyright by the 
Published by the 
of the Massachuse 
Printed in Cambrii 


[assachusetts Institute of Technology, 19' 2. 
larles Stark Draper Laboi atory, a division 
;s Institute of Technology 
tge, Massachusetts, U.S. v., 1972. 


ii 



PfiECEDINU hAo, 


^‘^ix A uj. 


A ILJlluL) 


E-2648 

CROSS FLOW IN A STARVED EHD CONTACT 


ABSTRACT 


It is suggested that fluid film thickness in a starved EHD 
contact is determined by the rate at which lubricant is forced from 
the Hertz zone normal to the rolling direction. A simple model 
allows calculation of this flow rate, together with related quantities 
of interest to the designer. 
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NOMENCLATURE 



Units 

Definition 

A, B, C, D 


constants 

E 

in. 

pitch diameter 

F 


tan' 1 ( 2H - K )] 

L (r - KH + H^) K y/5 J 

H 


h/h 0 

K 


(- Q- /Q a.) 1 ^ 3 
' y in v out' 

L 

in. 

half length of Hertz zone 

N 

lbf 

normal load carried by contact 

P 

lbf/in. 3 

mean Hertz pressure 

^out (in) 

in. 3 /s 

flow out of (into) Hertz volume 

T 

$ 

real time 

V 

. 3 
in. 

Hertz half volume 

Z 


number of balls 

d 

in. 

ball diameter 

f 


ratio between contact time and real time 

h 

in. 

film thickness 

t 

s 

contact time 

w 

in. 

Hertz width 

a 

in. 2 /lbf 

pressure viscosity coefficient 

V 

poise 

viscosity 

0 

rad 

contact angle 

Subscripts 



0 


refers to conditions at t s 0 

p 


refers to conditions prior to t * 0 

e 

1 

refers to conditions at t = co 
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CROSS FLOW IN A STARVED EHD CONTACT 


INTRODUCTION 


The fundamental elastohydrodynamic problem is to calculate lubricant film 
thickness in a contact where elastic deformations have the same order of magnitude as 
the film thickness. Various solutions are available, falling into two types: a 

Grubin or "inlet" analysis, which solves Reynolds equation in the inlet up to the 
edge of the Hertz zone (1,2, 3,4); and a "complete" analysis, involving simultaneous 
solution of the Reynolds and Hertz equations over the entire contact (5). In general 
any Grubin analysis ignores flows within the Hertz zone, and the available complete 
solutions are limited by mathematical complexity to a consideration of one-dimensional 
flow In the direction of rolling. Moreover, these solutions all suppose that film 
thickness is determined by parameters of load, speed, viscosity, geometry, and elas- 
ticity. The tacit assumption is made that thickness is not a function of the amount 
of lubricant, which is presumed to be present in excess at the contact. Attempts 
have been made to impose the effects of side leakage (e. g., 2) and oil supply limi- 
tation (e. g., 4) on the classical solutions. However, the modifications seem arti- 
ficial and do not change the underlying theoretical model. 

Considerable experience using ball bearings on the spin axis of precision gyro- 
scopes has shown thct successful operation can be obtained although the bearings are 
starved in an EHD sense. Starvation is defined as any operating condition such that 
an increase in oil available to the contact will result in an increase in film thick- 
ness (6). In this paper, starved EHD behavior is considered in several situations: 
when no new oil enters the contact, and when a step change in one of the rolling 
parameters is introduced at an equilibrium contact. Arguing that the important oil 
fiow within the Hertz zone Is across the rolling direction, one can formulate a simple 
theory which explains quantitatively some aspects of starved-ball -bearing behavior, 
e. g., the "oil jag." 
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OIL JAGS 


An oil jag in a gyro spin-axis bearing is detected as an abrupt increase in 
driving torque demand, followed by a decay lasting over many seconds back to the 
original torque level (7,8). If the gyroscope is being accelerated perpendicular to 
its spin axis during the jag (spin axis and output axis horizontal) , a simultaneous 
and one-one excursion, positive or negative depending upon which bearing of the pair 
is involved, is observed in its torque-to-balance signal (Fig. 1). This indicates an 
axial shift of the gyroscope CG. 


Driving 

Torque 


Torque 

to 

Balance 


t 



Fig. 1 . Disturbances in driving torque and torque to 
balance dur'ng oil jag. 
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Suppose, in one ball retainer, that oil flows in the centrifugal field to its 
OD, and accumulates in a droplet. The drop will grow to a cricital size determined 
by its surface tension and the field, whence it will fly off into the ball track. 
Energy required to roll the ball through the deepened oil is reflected in the driving 
torque excursion, and an axial shift of the gyro CG due to the increased film thick- 
ness in one bearing produces the torque-to-balance excursion (Fig. 2). 

Various other jag observations are exnlained by this model; stroboscopic 
observation has confirmed oil droplets on the retainer OD, together with a torque 
excursion upon their disappearance; artificial jags have been produced with a "jag- 
gun" which throws a known amount of oil into the ball track on tldmand. Thus the model 
appears correct, and associates a known increase in film thickness with a known 
increment in the amount of oil available to the contact. Further, it gives the time 
required to reduce the film back to its original thickness. 




increased h 
due to jag 



— axial CG 
displacement 


unbalance 

torque 



Spin Axis 


equilibrium h 
of steady state 


BEARING A 


BEARING B 


Fig. 2. Oil-jag-mechanism geometry. 
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FLOW MODEL FOR A STARVED CONTACT 


In a starved bearing, the oil on the track is about as thick outside the con- 
tact as in it. Stroboscopic observation of interference fringes confirms that films 
so thin flow very slowly in the centrifugal fields due to btaring rotation. The 
working oil is thus stationary wiui respect to the track outside the contact, and is 
unloaded most of the time. The oil momentarily in the contact may flow in the 
direction of rolling and also across the direction of rolling under the Hertz 
pressure acting on it. Its viscosity will be greatly increased while it is under 
this pressure. Flow in the rolling direction cannot alter the average film thickness 
since it removes no oi from the track. Only transverse oil flow can change the film 
thickness. The long-term decays in thickness observed with oil jags thus imply that 
transverse oil flow out of the Hertz zone should be considered in a starved EHD con- 
tact. 

To estimate this outflow, approximate the volume of oil in each half of the 
Hertz zone at any time as (Fig. 3) 

V = L w h [1] 

where L and w are the Hertz dimensions of the contact, and h is the film thickness. 
Calculate the oil flow out each end of this box due to a pressure drop equal to the 
mean Hertz pressure, P, of an oil of viscosity a /« a/(P) as 


Q = 


P h 3 w 

JTVT 


[ 2 ] 


Flow ir the roll direction may be neglected as it must average to zero. 

dV = L w dh - - Qdt [3] 




[ 4 ] 
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Fig. 3. Oil-film geometry in Hertz zone. 

inhere t is the tine an element in the race track is in the contact. It is a fraction 
f, of real tine T, given by the proportion between Zw, the total contact width, and 
fr{t + d cos B), the outer-race-track circ (inference. Thus, 

t = fT, f * Zw/(E + d cos d ) v [5] 

written for an outer-race contact in a bearing having Z balls of diameter d, pitch 
diameter E, and contact angle B. 

Finally, 



which relates the decay of a film of original thickness h Q to real time T in a bearing 
where no new oil enters the wear track. 

This calculation has been carried out for the R3-size bearing (Fig. 4) which 
produced the jags shown in Fig. 1. The results for two running conditions are plotted 
in Fig. S. It takes about three hours of real-time running to pump half the oil out 
of a contact initially 20 x 10" 6 inches thick under 150,000 lbf/1n. 2 . When the con- 
tact is 4 x 10~ 6 Inches thick, it takes about 100 hours to pump out half the remaining 




Fig. 4. R3 ball bearing. 

oil. These numbers are not out of line with life experienced in retainerless bearings 
which have no provision for the recirculation of oil. 

2 

If the decay is recalculated for a Hertz stress of 100,000 Ibf/in. , there is 
a decrease in the life of the contact. This is because of the very great decrease in 
the viscosity of the lubricant with pressure, from 10® to 10* poise in this case (9). 

THE STARVED EQUILIBRIUM CONTACT 

Since it is possible to run starved gyro bearings essentially indefinitely (10), 
tnere must be a flow of lubricant into the wear track to balance the outflow calcu- 
lated in the preceeding section. It has been confirmed experimentally (11) that oil 
circulates from a reservoir in the ball retainer to the track and back, although the 
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Fig. 5. Film decay without inflow. 


details of the process have not been explored. At equilibrium, therefore, there is 
an inflow equal to the transverse outflow, and an equilibrium film thickness. If 
equilibrium is disturbed by a change in some running condition, the transient response 
in film thickness at the contact may be calculated. 

In the following analysis, a step change is assumed to dccur at t * 0; sub- 
script p refers to equilibrium conditions prior to the change (t <0); subscript 
o to conditions at the change (t ■ 0); and subscript e to final equilibrium (t «®). 
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Thus, at an equilibrium contact for t < 0 


Q. = Q . 
y in p y nut p 


h w h ~ 

= tHjT' = Q in 


[7] 


P P 


It is supposed that tht inflow to the contact, Q in , is constant throughout the tran- 
sieit, since it is th»: result of long-term averaging processes in the retainer, but 
that the outflow ste^ at t 8 0 since it depends directly on contact conditions. For 
t >0 

S-3f[Wl‘>] =, in- ( W‘ M 




<W* * -f?Vr 0 


[9] 


so that 


12vA / — T s / dt 

. 12 ttAJL- - P_w h n 

h_ o o in oo 0 


[10] 


Equation [10] may be inte^ated to get 


where 


t w h n 
0 0 0 


15 


in 


K [F(H,K) - F(1,K). = t 


3 ^in 

K 3 = - , H = h/h 

g out_ 0 


[ 11 ] 


and 


F(H,K) 


[> in L.ya [12] 

\ £ (r - KH + tr) K/3 J 


F(H,K) has been plotted in Fig. 6 for relevant values of H and K. 
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F(H. K) 



H « h/h Q 


Fig. 6. F(H,K) versus H. 


FILM-THICKNESS STEP 

An oil jag nay be simulated as a step increase in film thickness from h p to h 0 , 
all other running parameters remaining constant. Then, from Eq. [12] 

vp- ■ H p “ | K | < 1 [13] 

Suppose Hp * - 0.5, meaning that the bearing jags to a new thickness twice the prior 
equilibrium value. Equation [11] becomes 

- 0.5(-S“^) [F(H, -0.5) - F(1 , -0.5)] - t [14] 

g in 
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Boundary conditions are satisfied since at t * 0 


F(H 0 , -0.5) = F(l, -0.5), H q * 1 and h s h Q [15] 

and at t = oo , F(H , -0.5) = -oo occurring (Fig. 6) when H = -K or 

C v 

0.5 = H a = h_/h , r\ a = 1/2 h„ [16] 

e e o’ e ' o 

Thus, for this case Eq. [11] describes a time decay from the stepped film thickness 
at t = 0 back to the prior equilibrium thickness at t = oo . 

As a numerical example, suppose an R3 bearing (F = 0.01) is running on an 
equilibrium film h p = 20 x 10”^ inches thick, with a Hertz pressure P p = 10^ lbf/in. 
and oil viscosity v = 1.1 x 10^ poise. At t = 0 the bearing jags to h Q = 24 x 10” 
inches, giving a K value of - = -0.833. Using Fig. 6 and Eq. [11], the decay 

shown in Fig. 7 is calculated. It may be compared with jags experienced in similar 
bearings (Fig. 1). 


Real Seconds 



Fig. 7. Film decay after jag. 




LOAD STEP 


Suppose that at t = 0 the load carried at a starved EHD contact steps from 
Np to N 0 . From elasticity theory 


(P, w, L) = (A, B, C)N 1/3 


[17] 


where A, B, and C are constants which may be calculated according to the particular 
bearing ur.der consideration. 

Assume an exponential dependence of viscosity on pressure, as is done in 
Ref. (4): 


* V 


*<v 


v 

H = u. 


A a(N 


1/3 


- »J' 3 > 


[18] 


Thus, at the load step, the Hertz pressure, the Hertz dimensions, and the viscosity 
all change to new values. However, the thickness does not step, hp = h Q . Its time 
dependence for t >0 is given by Eq. [11]. From Eq. [12] 


k 3 = - 


«1n. 

^out. 


P n wL n 
p p 0 0 

P w L n p 

0 0 p p 


[19] 


or, in terms of load 


N d 1/3 


A a (N 


1/3 


N ] / 3 ) 

P 7 


[20] 


The exponential part of K ensures that j K | >1, meaning that H values will be found 
on the left branch of the F curves in Fig. 6. To calculate the final equilibrium 
fPta thickness, note that at t ■ o° , -F(H,K) * oo. From Fig. 6 


■K = H 


e* 


[ 21 ] 


h fi may also be calculated by observing that at t =oo , Q Qut = Q^. Substitution 

‘'e 

gives the same result for h g . 

For this case, Eq. [11] thus describes a build-up In film thickness with time 
after a step Increase in load (Fig. 8). 


11 




Fig. 8. Film thickness history after load step. 


DISCUSSION 


This analysis indicates that to calculate the EHD film thickness in a starved 
contact, it is necessary to know the rate at which lubricant is supplied to the con- 
tact. In a conventional instrument bearing, oil is stored in the ball retainer and 
metered out Dy unknown processes at unknown rates at the different ball pockets. The 
direct calculation of film thickness in these bearings is therefore not feasible. 

Testing of retainerless bearings which incorporate a rational lubrication 
supply system has been started. The oil supplies for these bearings have been speci- 
fically designed to provide well-defined inflows. As such information becomes avail- 
able, thickness predictions from the cross-flow analysis can be tested. 

The model for cross flow is crude but is believed to reflect some reality. 
Calculation of life and jag decay seems adequate, and special experiments to check 
some other predictions are under consideration. Predictions include independence of 
jag decay period with speed, and increase in decay period with stress. Depending on 
the outcome of such experiments, this model can be improved or scrapped. 
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